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ABSTRACT 

We present a census of blue horizontal branch (BHB) and blue straggler (BS) stars be¬ 
longing to dwarf galaxies and globular clusters, and compare these counts to that of the Milky 
Way stellar halo. We find, in agreement with earlier studies, that the ratio of BS-to-BHB stars 
in these satellite populations is dependent on stellar mass. Dwarf galaxies show an increasing 
BS-to-BHB ratio with luminosity. In contrast, globular clusters display the reverse trend, with 
Wbs/.^bhb 1) decreasing with luminosity. The faintest (L < lO^L©) dwarfs have similar 
numbers of BS and BHB stars (iVBs/.^BHB ~ !)> whereas more massive dwarfs tend to be 
dominated by BS stars {Nbs/Nbhb '^2 — 40). We find that the BS-to-BHB ratio in the stellar 
halo is relatively high (Wbs/^bhb ~ 5 — 6), and thus inconsistent with the low ratios found 
in both ultra-faint dwarfs and globular clusters. Our results favour more massive dwarfs as 
the dominant “building blocks” of the stellar halo, in good agreement with current predictions 
from ACDM models. 
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1 INTRODUCTION 

The Milky Way is a cannibal; throughout its lifetime it captures and 
destroys smaller dwarf galaxies. The remains of destroyed dwarfs 
are splayed out in a diffuse stellar halo, while the dwarfs evading 
destruction comprise the satellite population that orbits the Galaxy. 
Despite this well-established, generic picture of stellar halo forma¬ 
tion, we have very little understanding of what the building blocks 
of the halo actually are; is the halo built up from many small mass 
tidbits, or from one (or two) massive dwarf(s)? 

The chemical properties of halo stars have often been used 
to connect them to their progenitor galaxies. For example, the re¬ 
lation between [a/Fe] and [Fe/H] is an indicator of the rate of 
self-enrichment, and therefore can be linked to the host galaxy’s 
mass. However, the [a/Fe] abundances of halo stars appear to dif¬ 
fer significantly fr om those of the (classical) dwarf galax y satellites 
in the Milky Way jTolstov et alJl200^ ; IVenn et alj|2004h . whereby 
the halo stars are typically more a-enhanced at a given metallicity. 
Thus, there is little evidence for the accretion of fragments similar 
to the present-day dwarf spheroidal population. 

The mismatch in chemical properties between the bulk of the 
halo stars and the stars belonging to dwarf spheroidals can perhaps 
be reconciled if the Milky Way halo progenitors are biased towards 
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massi ve, early accretion events ^Robertson et aljBoOSi ; iFont et all 
1 20061) . The combination of high-mass and early accretion, can lead 
to abundance patterns (at least in the [a/Fej-lFe/H] plane) similar to 
that exhibited by the present day halo stars. This scenario has been 
supported by recent evid ence of a “break” in the stellar halo den sity 
profile at r ^ 25 kpc jPeason et al. l l201ll;ISesar et alj|201 ll) . In 
IPeason et alj l l2013l) . we argue that this break could be evidence 
for a major (relatively early) accretion event. However, this is not a 
unique solution; the same broken profile can plausibly be produced 
from multiple, but synchronized, lower-mass accretion events. 

A different scenario posits that analogues of the “ultra-faint” 
dwarf galaxies could contribute significantly (a t least at the metal- 
poor end) to the present-day st ellar halo (e .g. Frebel et al. IH; 
IClementinill2010l) . For ex ample, Clementinil 
Oosterhoff classification ( IOosterhofllll93S 
ultra-faint dwarfs is in better agreement with the stellar halo com¬ 
pared to the more massive dwarffl Thus, an alternative view is that 
the stellar halo is built-up from a very large number of puny dwarfs. 
Finally, bear in mind that the characteristic building blocks of the 
stellar halo need not be dwarf galaxies. Previous work has argued 


(e.gjFreb 

oil yOiy) argue that the 
of RR Lyrae stars in 


^ However, we note that iFiorentino et alj j20l4) recently showed that 
the period and luminosity amplitudes of RR Lyrae stars in the halo ai'e 
more consistent with massive dwarfs (such as Sagittarius) than lower-mass 
dwarfs. 
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that a significant fraction of the stellar halo ( up to 50%) could be 
assembled from de stroyed globular clusters dCarretta et al.ll20f^ : 
iMartell etalJl201lh . 

Despite the wealth of work attempting to decipher the mass 
spectrum of accreted substructures, we currently lack a clear pic¬ 
ture of what made up the stellar halo, and when. In this letter, 
we use an alternative approach to gain insight into the progeni- 
tors of the Galactic h alo. Recent work bv iMomanvl l l2014h (see also 
iMomanv et al.l200% showed that the number ratio of blue straggler 
(BS) to horizontal branch (HB) stars in dwarfs and globular clusters 
is dependent on the satellite’s stellar mass. Hence, this ratio could 
potentially be used to constrain the mass spectrum of substructures 
that contributed to the stellar halo. With this aim in mind, we pro¬ 
vide a careful comparison between the number ratio of BS-to-blue 
horizontal branch (BHB) stars in different Milky Way companions 
(classical dwarfs, ultra-faint dwarfs and globular clusters) and the 
stellar halo overall. 


2 A-TYPE STAR POPULATIONS IN THE MILKY WAY 
HALO 

In this section, we identify the BHB and BS pop ulations in dwarf 
gala xies, globular cluster s and the stellar halo. iMoman 3 l20li) 
(also lMomanv et al]|2007h showed that the ratio of BS to HB stars 
varies as a function of luminosity for satellites in the Milky Way. 
However, in their study the entire HB was considered, which in¬ 
cludes the red horizontal branch (RHB) and the extended blue tail 
of the HB. The RHB is notoriously difficult to identify in the stel¬ 
lar halo, and current BS-to-HB ratios in the stellar halo are upper 
limits as only BHB stars are included. Hence, in this work we con¬ 
sider the BS-to-BHB ratio for a fair comparison between satellites 
and the field halo. Our use of BHB stars on the denominator of this 
population ratio could be perceived as problematic, particularly if 
the BHB population is scarce, or does not exist at all in some satel¬ 
lites. However, it is worth pointing out that, to our knowledge, there 
isn’t a single dwarf galaxy that does not have any BHB stars. On 
the other hand, some very metal-rich globular clusters are devoid 
of a BHB population, and we discuss this further in Section l2^ 
We note that our choice of BS-to-BHB ratio as a probe of 
the stellar halo progenitors is made for both physical and practi¬ 
cal reasons. The BS-to-BHB ratio is arguably the cleanest popu¬ 
lation relation that can be measured in both satellite galaxies and 
the stellar halo (see Section l23t . In particular, redder populations 
such as RHB and red giant branch (RGB) stars suffer from severe 
foreground contamination, and are much more difficult to isolate in 
the halo with photometry alone. However, the main advantage of 
using these A-type star populations is that the BS-to-BHB ratio is 
easier to quantify in the stellar halo than the total numb er of BHB, 
BS, RGB, RHB etc. stars alone (see lPeason et al.l201 ll and Section 
[23] ). 

2.1 Dwarf Galaxies 

Our compilation of dwarf galaxies in the Milky Way is obtained 
from a variety of photometric data sources in the literature (see 
Table [T]). We ensure that our sample only includes photometric 
data deep enough to reliably identify the BS population (typically 
~2 magnitudes fainter than BHBs) from the colour-magnitude di¬ 
agram (CMD), and we only include datasets where Wbhb > 1 and 
(Vbs > 1- This excludes some of the more distant dwarfs without 
sufficiently deep photometry (e.g. Canes Venatici I), and some of 
the ultra-faint dwarfs with very few stars (e.g. Segue I). 


Our dwarf sample excludes cases with kn own recent star for - 
mation (e.g. Fornax, Leo I, Carina - see e.g IWeisz et alJl^Oldl) . 
where contamination by young stars inhibits reliable estimates 
of t he BS population. Very young stars (~ 1 — 3 Gyr, see 
e.g. ISantana et al.l l2013l) can mimic BS stars in dwarf galax- 
ie s, and we a r e guid ed by the star formation histories derived 
in IWeisz et^ ) |2014|) to exclude these cases where possible. For 
consistency, we convert all magnitudes into SDSS bandpasses. 
Johnson-Cousins m agnitudes ar e conve rted to gri SDSS filters us¬ 
ing the relations in jjordi et ^ )l2006l) . and HST/ACS filters are 
converted i nto Johnson-Cousins bandpasses using the procedure 
outlined in ISirianni et al.l l l2005l) . The magnitudes and colours we 
use have been correc ted for extinction following the prescription of 
ISchlegeletaLl l ll998h . 

For a fair comparison with the stellar halo (see below), only 
A-type stars with —0.25 < g — r < 0 are used. In cases where g — i 
colour is most appropriate (e.g. for V, I filters), we used bright A- 
type stars from SDSS (16 < g < 17) to calibrate a linear relation 
between g — i and g — r. We find that the colour range —0.25 < 
g — r < 0.0 roughly corresponds to —0.44 ^ g — i ^ —0.11 for 
A-type stars. We use the globular cluster sample (see below) with 
gri photometry to ensure that our selection of A-type stars in g — i 
is consistent with our selection using g — r. 

Some example CMDs are shown in Fig. [T] The selection re¬ 
gion for BHB and BS stars are shown with the blue and red poly¬ 
gons, respectively. We use the Trilegal Galaxy model jGirardi et al.l 
1 20051) to estimate the foreground contamination included in our A- 
type star samples. The estimated foreground in the BHB and BS 
CMD selection regions is subtracted before the BS-to-BHB frac¬ 
tions are computed. Note that in some cases control-fields are avail¬ 
able, and we use these to ensure that our estimated contamination 
from the Trilegal model is doing a reasonable job. In general, the 
contamination in the blue (g — r<0) region of CMD space probed 
in this work is minimal. 

The resulting BS-to-BHB ratios are given in Table [T] and 
shown in Fig. The quoted error estimates only include Pois¬ 
son noise. For datasets where we are privy to the full photometric 
error distribution, we find that our measurements are not signifi¬ 
cantly affected by photometric uncertainties, and in most cases, the 
error budget is indeed dominated by number statistics. There are 
several unavoidable sources of error apparent when computing the 
BS-to-BHB ratio: (i) uncertain foreground/background subtraction; 
(ii) confusion between BS stars and normal main sequence stars; 
and, (iii) radial gradients in dwarfs. However, the general agree¬ 
ment between BS-to-BHB ratios from different data sources (dif¬ 
ferent FOV, filters, sample size etc.) of the same dwarf is encour¬ 
aging, and suggests that these potential systematic uncertainties are 
not significantly affecting our results. Where there are multiple data 
sources for the same dwarf we show the weighted (by inverse vari¬ 
ance) mean value of Abs/Abhb in Fig.|^ 

2.2 Globular Clusters 

We also show in Fig. the BS-t o-BHB number ratio for globu¬ 
lar clusters in the [An et'H] l l2008l) sample. These globular clusters 
have SDSS photometry and the BHB and BS populations are iden¬ 
tified from the CMDs in the same way as the dwarf galaxies. We 
only include globular clusters with Abhb > 1 and Abs > 1, and 
ensure that the photometry is deep enough to identify the BS popu¬ 
lation. This leaves a sample of 12 globular clusters that satisfy our 
requirements. An example of a globular cluster CMD is shown in 
the right-hand panel of Fig.[T] 
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Figur e 1. Three example colour-magnitude diagrams in gri SDSS filters, with original photometry from iBelokurov et alj j20o3)jHel j j2Q05h and I An et alj 
<200 SI) , respectively. The selection of BHB/BS stars are indicated with the blue/red lines, respectively. 


Name 

Mv 

Photometry 

FOV 

Nbs / Nbhb 

Ref 

( Nbs / Nbhb ) 

Bootes I 

-6.3 

Blanco/Mosaic-II (g, i) 

36' X 36' 

1.3 ±0.4 

B06 

1.5 ±0.5 



Subaru/Suprime-Cam {V, I) 

34' X 27' 

2.5 ±0.9 

012 




HST/ACS (F606W, F814W) 

(5) 3.4' 2 

3.0 ±2.5 

W15 


Canes Venatici II 

-4.9 

Subaru/Suprime-Cam (g', i') 

34' X 27' 

0.6 ±0.3 

B07 

0.6 ± 0.2 



HST/ACS (F606W, F814W) 

3.4' 2 

0.7 ±0.4 

W15 




HST/WFPC2 (F606W, F814W) 

2.4' 2 

0.7 ±0.6 

H06 


Cetus 

- 11.2 

HST/ACS (F475W, F814W) 

3.4' 2 

45.4 ±11.5 

M12 

45.4 ±11.5 

Coma Berenices 

-4.1 

Subaru/Suprime-Cam (g^ i^) 

34' X 27' 

0.5 ±0.4 

B07 

0.5 ±0.4 

Draco 

-8.8 

INT/WFC (V, I) 

~ 1 deg2 

5.4 ±1.2 

AOl 

5.3 ±1.1 



HST/ACS (F555W, F814W) 

3.4' 2 

4.0 ±3.2 

W15 


Hercules 

- 6.6 

HST/ACS (F606W, F814W) 

3.4' 2 

1.0 ± 0.6 

W15 

1.0 ± 0.6 

Leo II 

-9.8 

HST/WFPC2 (F555W, F814W) 

2.4' 2 

9.2 ±2.7 

H05 

10.0 ± 1.9 



HST/ACS (F555W, F814W) 

3.4' 2 

14.8 ±4.9 

W15 




HST/WFPC2 (F606W, F814W) 

2.4' 2 

9.4 ±2.7 

H06 


Leo IV 

-5.8 

Subaru/Suprime-Cam {V, I) 

34' X 27' 

3.2 ±1.5 

012 

1.7± 1.0 



HST/ACS (F606W, F814W) 

3.4' 2 

1.0 ± 1.0 

W15 


Sagittarius 

-13.5 

MPl/WFl (V, 1) 

~ 1 deg^ 

10.0 ± 1.0 

M03 

10.0 ± 1.0 

Sculptor 

- 11.1 

MPI/WFI (B, V, I) 

34' X 33' 

2.6 ± 0.1 

R03 

2.0 ±0.4 


CTIO/MOSAIC (V, I) 

4 deg^ 

1.8 ± 0.1 

Til 


Sextans 

-9.3 

CFHT/CFH12K (B, V, I) 

42' X 28' 

6.2 ± 1.1 

L03 

6.2 ± 1.1 

Tucana 

-9.5 

HST/ACS (F475W, F814W) 

3.4' 2 

3.7 ±0.3 

M12 

3.7 ±0.3 

Ursa Major I 

-5.5 

Subaru/Suprime-Cam (V, I) 

34' X 27' 

0.3 ±0.2 

012 

0.4 ±0.2 



INT/WFC (B, r) 

23' X 12' 

1.5 ±1.0 

W05 


Ursa Minor 

-8.8 

INT/WFC (B, R) 

0.75 deg2 
( 2 ) 2 . 4 ' 2 

1.0 ± 0.1 

C02 

1.0 ± 0.1 



HST/WFPC2 (F555W, F606W, F814W) 

2.3 ±1.0 

H06 



Table 1. The dwaif galaxies used in this work. We list the dwaif name, absolute visual magnitude, the photometry used to calculate population ratios, 
approximate FOV, B S-to-B HB number ratio , appro priate references to the photo metric data sources, and (w eighte d) average B S-to- BHB number ratio. API : 
lAparicio et al.l 

' '^ donaco et alJ<20030 . 012: IUkamoto et al.H2U12l) . K.03!lRizzi t 


L03: lLee et SLooT MO^ iM 
wn.^i lwillmaii e.t 311120051) 


, _ . _; photo m( _ _ ^ _ 

200 1[) . B06: lBelokurov et alH2006h . BO VilBelokurov et al.1 J2007h . C 02:ICan-era et alj|2002l) HQ5:lHeldl2005h. H 06: lHolt7man e.t al] j2006h . 

^^~~ ~eraiTlM03h . O12: lokamoto‘etaiT72012h . R03: lRizzi eTalTl^OO^ ), T1 l: lde Boer et alJl201 ih . W15: Weisz et al. in prep.. 


We note that we do not include the relatively metal-rich 
([Fe/H] > —0.8) globular cluste rs that do not have a BHB popula¬ 
tion, but do have BS stars (see e.g IPiotto et al.l2003l . These systems 
would boast abnormally high BS-to-BHB ratios and could poten¬ 
tially contribute BS stars to the halo. However, gi ven the metal-poor 
nature of the stellar halo (([Fe/H]) ~ —1.5. Urozic et ahlboOSl : 
lAnetalJ2013h it is reasonable to assume that these metal-rich glob¬ 
ular clusters are not significant contributors. 


2.3 Stellar Halo 

The identification of BS and BHB stars in the stellar halo is not 
as straightforward. While at bright magnitudes {g < 18.5), A-type 
stars can easily be distinguished from white dwarfs and quasars us¬ 
ing ugr photometry, BS and BHB stars cannot be cleanly separated 
using photometry alone. 

In IPeason et al.l ( 1201 ih (hereafter, DBEll), we used A-type 
stars selected from SDSS to measure the density profile of the stel¬ 
lar halo out to D ~ 40 kpc. DBEl 1 took advantage of the overlap- 
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Figure 2. Number ratio of BS-to-BHB stars as a function of absolute vi¬ 
sual magnitude. Dwarfs and globular clusters are shown with black squares 
and green circles, respectively. The quoted error bars only include Poisson 
noise. The halo number ratio is computed using the results of DBEl 1; the 
derived halo density profile is used to convert the overall number ratio (in 
a fixed magnitude slice) to a number ratio at fixed volume. The solid gray 
and line-filled green bands illustrate the correlation between AIbs/^bhb 
and Mv for dwarfs and globular clusters, respectively. 

ping, but distinct, ugr distributions of BS and BHB stars (see Fig. 
2 in DBEll). The BFIB and BS populations were modeled simul¬ 
taneously with class probabilities based on ugr photometry alone. 
This method resulted in two quantities important for this work: 1) 
an estimate of the number ratio of BS-to-BffB stars in a fixed mag¬ 
nitude slice (see Table 1 in DBEl 1) and 2) a measure of the stellar 
halo density profile, under the assumption that both BHB and BS 
populations follow the same density profile. 

In order to compare the stellar halo with the satellite popula¬ 
tions, we must take into account the different volumes probed by 
BHB and BS stars in a fixed magnitude slice (BS stars are ~ 2 
mag fainter than BHB stars). Thus, we use the ratio Pbs/Pbhb> 
where p^g = Nbs/Vbs and p^ns = -^bhb/Vbhb- Here, A^bs 
and A'bhb are the numbers of BS and BHB stars in a fixed mag¬ 
nitude slice, and the volumes (Vbs, Vbhb) are given by equation 
(9) in DBEll. The resulting ratios are 4.9 ± 0.1 and 6.4 ± 0.1 
when stars belonging to the Sagittarius stream are excludecQ or in¬ 
cluded, respectively. The error estimates take into account the dif¬ 
ferent likelihoods of stellar halo density models. The halo ratios are 
shown with the purple and blue lines in Fig.|2 


3 POPULATION RATIOS: COMPARING SATELLITES 
WITH THE STELLAR HALO 

Our compilation of BS-to-BHB number ratios for halo populations 
is shown in Pig. as a function of absolute magnitude. Dwarfs 
and globular clusters are displayed with the solid black squares and 

^ Using the same mask defined in DBEl 1. 


open green circles, respectively. The BS-to-BHB number ratios in¬ 
crease with luminosity for the dwarfs, but the opposite trend is seen 
for t he globular clust ers. 

|m Oman ^( l20li) showed that the BS-to-HB ratio for dwarfs 
decreases with absolute magnitude. We find the opposite trend for 
dwarf galaxies when only the BHB pop ulation is includ ed on the 
denominator. This difference is because iMomanvl ( l2014h includes 
all HB stars (BHB, RHB and the extended blue tail of the HB) 
in their analysis, so their trend is likely due to the more massive 
dwarfs having a more prominent RHB. As stated earlier, the RHB 
population is extremely difficult to quantify in the stellar halo, so 
the BS-to-BHB ratio provides a more robust comparison between 
satellites and halo stars. 

The difference in trends shown for GCs and dwarf galaxies 
is likely related to the different BS formation mechanisms in these 
systems. The two m ain established routes of BS production (see e.g 
IPavies et al.l2004h . from collisional-binaries and primordial, wide- 
binaries, have different significances in dwarfs and clusters; both 
formation channels act in globular clusters, whereas the low stellar 
density environments of dwarf galaxies precludes the occurrence 
of collisional binaries. Additionally, the higher densities (and col- 
lisional probabilities) in more massive globular clusters can lead to 
BS disruption (through 3-body interactions), but this process is not 
important for the (similar mass) dwarfs. This likely explains the 
large differences in Nbs / Nbbb fractions at Mv ~ — 9 between 
dwarfs and globular clusters. 

The fainter dwarfs {Mv > —7.5) have similar BS-to-BHB 
number ratios to globular clusters at comparable luminosities, 
whereas more massive dwarfs have much higher ratios than globu¬ 
lar clusters. 

While the brighter Milky Way dwarfs have much larger BS- 
to-BHB number ratios than the fainter dwarfs, there is also a good 
deal of scatter. Por example. Sculptor and Cetus have very simi¬ 
lar absolute magnitudes {Mv ~ —11) but very different number 
ratios, Nbb/Nbbb ~ 2 for Sculptor and Nbb/Nbbb ~ 40 for 
Ce tufl The sta r form ation histories of these two dwarfs derived 
by IWeisz et"^ ( l2014h from HST photometry are also very differ¬ 
ent, where Sculptor has a much older stellar population. It is clear 
that at fixed luminosity the star formation histories (and hence BS- 
to-BHB number ratios) can vary substantially, especially for more 
massive dwarfs. 

Despite the large scatter for bright dwarfs, it is clear 
that the low BS-to-BHB number ratios for ultra-faint dwarfs 
{Nbs/Nbbb ~ 1) and globular clusters (Wbs/Abhb < 1) are not 
compatible with the relatively high BS-to-BHB ratio in the Milky 
Way stellar halo. Thus, it is unlikely that the bulk of the stellar 
halo was built up from (a very large number of) low luminosity 
systems such as ultra-faint dwarfs and/or globular clusters. This is 
in agreement with the current model predictions from ACDM sim¬ 
ulations, postulating that stellar halos are generally dominated by 
massive accretion events jBullock & Johnst^l2005l : ICooner et al.l 
l20ld ; lDeason et alj2013l) . We do note, however, that if the progen¬ 
itor satellites were drastically different to the surviving populations 
today, then we must be more circumspect regarding our compar¬ 
ison with halo stars. For example, globular clusters destroyed a 
long time ago (~ 10 Gyr) may not have had time for collisional 
processes to occur, and thus the BS population may be very dif- 


® The unusually high BS-to-BHB number ratio in Cetus may result from 
contamination by young (1 — 2 Gyr) stars. Yet, to our knowledge, there is 
no evidence for such population in the dwarf. 
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ferent in these pro to-clusters. On the other hand, recent work by 
iBrown et alj ( l2Q12h arguing that the UF-dwarfs are predominantly 
ancient (~ 12 — 14 Gyr) populations, suggests that we are not sig¬ 
nificantly biased when comparing with the “survivors” at these low 
mass-scales. 

4 CONCLUSIONS 

In this Letter, we compiled a sample of BS and BHB stars in dwarf 
galaxies, globular clusters and the Milky Way stellar halo with the 
aim of comparing the BS-to-BHB ratio for different halo popu¬ 
lations. We ensure that our selection of BS and BHB stars is as 
consistent as possible (i.e. using the same photometric system and 
colour cuts) between different datasets, and correct the approxi¬ 
mate number ratio of BS-to-BHB stars in the stellar halo (at fixed 
magnitude slice) for volume effects. Our main conclusions are as 
follows: 

• The number ratio of BS-to-BHB stars in dwarf galax¬ 
ies increases with increasing luminosity. Ultra-faint dwarfs have 
Nbs/Nbhb ~ 1, while more massive dwarfs can range from 
Nbs /Nbhb ~ 2 to Nbs /Nbhb ~ 40. The large scatter for more 
massive dwarfs is probably due to the wide variation in star forma¬ 
tion histories. 

• GCs tend to have low BS-to-BHB ratios, Nbs/Nbhb (S 
which decreases with increasing luminosity. The different trends 
shown by globular clusters and dwarfs likely reflect the different 
for mation mechanisms of BS stars in these two populations (see 
e.g. lSantana et al .1120131: lMoman^l2014h . 

• The relatively high BS-to-BHB ratio in the stellar halo 
(Nbs /(Vbhb ~ 5 — 6) is inconsistent with the low ratios found for 
ultra-faint dwarfs and globular clusters. This result argues against 
ultra-faints and globular clusters being the dominant “building 
blocks” of the stellar halo, and instead favours more massive dwarfs 
as the more predominant progenitors. 
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